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We study the decay process in an open system, emphasizing on the relevance of the environment's 
spectral structure. Non-Marlcovian efltects are included to quantitatively analyze the degradation 
rate of the coherent evolution. The way in which a two level system is coupled to different envi- 
ronments is specifically addressed: multiple connections to a single bath (public environment) or 
single connections to multiple baths (private environments). We numerically evaluate the decay 
rate of a local excitation by using the Survival Probability and the Loschmidt Echo. These rates 
are compared to analytical results obtained from the standard Fermi Golden Rule (PGR) in Wide 
Band Approximation, and a Self-Consistent evaluation that accounts for the bath's memory in cases 
where an exact analytical solution is possible. We observe that the correlations appearing in a public 
bath introduce further deviations from the PGR as compared with a private bath. 
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I. INTRODUCTION 

Quantum Information Processing (QIP) requires an efficient and precise control on quantum dynamics jll]- Deco- 
herence conspires against this objective as it leads to progressive and smooth destruction of the quantum interferences 

t within a characteristic time t^. The main source of decoherence for solid-state spin systems such as quantum dots 
^ 6], donors in silicon [7|, defects in diamond [8l-[To|. and solid state nuclear magnetic resonance (NMR) [TT| - [l3| is 
the uncontrolled spin bath environment. While one can attempt different strategies such as quantum error correction 
protocols [13, [3l and dynamical decoupling [iB] , their specific efficiency depends on a deep understanding of how 
the environment behaves. Thus, understanding and mitigation of decoherence is one of the current challenges for 
quantum science and technology [l7l - [T9| . 

The interaction rate \/ts-£ between a quantum system S and its environment £ is typically evaluated from a 
Fermi Golden Rule (PGR), which assumes that the environment has a Markovian nature. Also, the interplay between 
the system time scales (e.g. and that of the interaction {ts-s) could result in striking effects. While weak 

interactions (\/ts-e ^ 2a;o) simply degrade dynamical interferences at a rate l/r^ oc l/r^-g, stronger ones may 
change the system's response radically, leading to a quantum dynamical phase transition in its dependence on 1/ts-s 
[l3| . Indeed, the possibility of a non- analytic behavior, e.g. in the energy spectrum [2Ci], appears because the system's 
effective Hamiltonian is non-Hermitian [2lll22l|. This, in turn, can be traced back to the fact that the environment has 
a number of degrees of freedom N which can be considered infinite. As P.W. Anderson put forth: "more is different" , 
and new physical phenomena may appear when this thermodynamic limit (N — > oo) is properly taken (23j . Indeed, 
there are several models of both non-Markovian and Markovian environments [12, [23 - [26j which show this dynamical 
phase transition as a function of the interaction strength. 

While the Markovian approximation is sufficient for most traditional applications, it leaves aside important memory 
effects and interferences in the time domain. These result from the coherent interaction between iS and f , and are 
becoming a topic of increasing interest The system-environment dynamics may go through different temporal 

regimes: a quadratic decay at very short times yields to the usual exponential PGR decay and much later it is 
transformed into an inverse power law decay (see Ref. [28|). Here, we will focus on the more relevant exponential 
regime, treated both with the PGR and with a Self- Consistent (SC) PGR which accounts for bath memory effects. 
Once again, a deep understanding of the environment's dynamics is central to identify the different regimes and to 
foresee possible dynamical transitions. 

A natural way to quantify the decoherence time is through the degradation of specific interferences, e.g. Rabi 
oscillations 29, 30] or mesoscopic echoes [3ll - l33 |. Alternatively, the implementation of a time reversal procedure, the 
Loschmidt Echo (LE) , allows the evaluation of the decoherence time by measuring the reversibility of the system's 
dynamics in presence of an uncontrolled environment. The LE can be accessed experimentally in spin systems [SGllSTj. 
confined atoms [sl], microwave excitations [s^, etc., and has become a powerful tool for quantifying decoherence, 
stability and complexity in dynamical processes in several physical situations [i^, |4l| . 

In this article we have the aim to quantify the role of bath's memory as well as specific correlations in the S ~ £ 
interaction on decoherence. With this purpose, we consider the evolution of two coupled spins in the fiand |t 
configurations in the presence of different spin environments with a fully characterized coherent dynamics [42|. The 
system's Rabi oscillations |43|] can act as a SWAP gate, and after appropriate mappings, this boils down to an 
excitation that jumps between two degenerate states A and B with a rate determined by the coupling "constant" 
Vab, which can be switched at will. Starting on state A, the return probability oscillates with wq = 2VyiB/?i, the Rabi 
frequency. In order to understand the incidence of environmental memory effects in the rate l/r^, we must consider 
the relation between the system's time scale (typically ruled by K/Vab) and the bath's inner excitation spreading time 
scale as determined by the density of directly connected states hNi , which, in general, is just a Local Density of States 
(LDoS) [44]. Also important is the specific form in which the system couples to the bath: each site may be coupled 
to a different environment (private bath) or both sites could share the same environment (public bath), and this 
would allow different correlations. A quantitative comparison between rates in these cases will enable a qualitative 
interpretation in terms of the bath's spectral structure and the effects of public and private S — £ interactions. This 
should deepen our understanding about how the involved correlations modify the degradation rate l/r^. Indeed, if 
the system's and the bath's time scales are similar, the problem cannot be treated within the Markovian paradigm 
of a "slow variable interacting with a fast equilibrating background" . It requires to be carefully addressed beyond 
the PGR. The appearance of mechanisms of correlated coupling in public environments has been previously pointed 
in the literature of open quantum systems, particularly in terms of the spin-boson models 45,, .46,] . It was shown 
(46j , at least for a simple model of a N— qubit register, that the decoherence increases linearly with N for independent 
reservoirs, while it grows with the square N'^ for a collective environment. This suggests that such a public interaction 
may lead to a strong amplification of decoherence. On the other hand, it has been pointed in the literature of 
error correction protocols that symmetric (public) S — £ coupling can be exploited to design states that are hardly 
corrupted by such a coherent environmental noise [47]. Also, in recent years the role of a common environment in the 
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entanglement correlations within a system [48|, |49| , and the possibility of creating and manipulating those correlations 
by environment-mediated interactions (soj have been explored. 

It is important to notice that spectral correlations within the bath, as described above, may become quite cumber- 
some to describe in statistical terms. However, realistic Hamiltonian models of the bath, as proposed here, should allow 
a simple and natural description of such correlations. Thus, the cases we analyze here can be casted directly into 1-D 
spin systems interacting by means of a planar (XY or flip-flop) or as chains of spins interacting with a double-quantum 
(flip-flip / flop-flop) effective Hamiltonians. Indeed, the spin-fermion mapping provided by the Jordan- Wigner Trans- 
formation (JWT) ^Sli] has been successfully exploited to p redict spin polarization dynamics in linear chains and rings 
[33} and results in full agreement with experiments (sj jslj. Furthermore, it allowed addressing polarization spin 
dynamics in homogeneous chains [2^ HO, [11] and dynamics of multiple spin coherences (ssI - Isg} by mapping them to 
fermionic excitations. In both cases, the local excitation is identified with a single fermion propagating in a tight- 
binding linear chain, and one can assume that the environment is described by an identical chain. Additionally, the 
dynamics of a tight-binding model can be mapped, sometimes quite straightforwardly, to describe the propagation of 
excitations in several scenarios like classical and quantum coupled oscillators plasmonic wave guides 
propagation [59'], spin-boson interactions and other models used for decoherence in quantum information 
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A fundamental issue about integrability of the chosen models is that they enable analytical expressions valid in the 
thermodynamic limit of an infinite number of spins. In fact, an analytically tractable bath means the capability to 
sum up an otherwise divergent perturbation series into a complex self-energy through the Dyson Equation, which only 
then is considered in the FGR approximation. Additionally, having a smooth LDoS, which arises from a continuous 
spectrum (N — > 00), avoids spurious resonances that might appear in finite systems. Besides, we will see that a well 
defined curvature of the LDoS will be central to quantify the corrections to the FGR. 

This paper is organized as follows. In Sec. |n]we present the S — £ tight-binding models that yields a fully solvable 
quantum dynamics (the underlying spin-fermion mapping is summarized in the Appendix In Sec. IIIIl the notions 
of Survival Probability and local Loschmidt Echo are presented on the face of their numerical implementation. We 
deal analytically with every case by working with the exact Green's Function in each case. Further details related to 
the GF's poles are discussed in Appendix [Bl In Sec. IIVI the analytical decay rates obtained with the FGR and the 
SC-FGR are presented for each case under consideration. These rates are compared with those computed numerically 
from the dynamics, contrasting the results of the Survival Probability (SP) and the Loschmidt Echo (LE) of the local 
excitation. All these magnitudes are analyzed under the light of the knowledge of the environment memory effects 
as provided by its spectral structure (LDoS). In the last section, further discussions and conclusions are presented 
to argue how the public bath is more effective to distort the decay process away from the usual Fermi Golden Rule 
(energy independent rate). 



II. TIGHT BINDING MODEL FOR EXCITATION DYNAMICS 



As pointed above, while our motivation lies mainly on spin dynamics under flip-flop interactions, we use the spin- 
fermion mapping to cast it in terms of tight-binding models that apply to a wide variety of systems. Thus, we leave 
to Appendix |X] a brief outline of how this mapping is achieved. Here, we present the models and analyze them in 
terms of straightforward single particle physics. These are basically variations of tight-binding infinite linear chains 
where the bath's memory can be fully characterized. Of course, the cases where the interactions network topology 
has branching points or loops would preclude the simple back transformation into spin systems. However, even in 
these situations some of the physics of the memory would remain. The general situations are sketched in FiglTJ 

For the cases considered in this article, the whole Hamiltonian can be written as: 

H = Hs+J2^£.^+^s-s, (1) 

where: 

Hs = Eac^a^a + Ebc'bCb - Vab{c\cb + c^c^), (2) 

with c\ and Cg (s € {A, _B}), the creation and destruction operators for fermions and Vab is the hopping amplitude 
that defines the system's only dynamical time scale provided that Ea = Eb- Each bath £1, (i/=l,2), with a spectral 
bandwidth of 41/, is described by the Hamiltonian: 

00 
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FIG. 1. (Color Online) Non interacting fermion modelization. I) A single site state connected to a semi infinite chain. II) One 
site coupled to an infinite chain. Since Vo <^V, both cases can be treated by a simple decay process, given by the FGR. Ill) 
Two sites; the initial configuration is given by a particle in the site A. The sites are coupled to a semi infinite chain by only 
one of them. IV) Same as III), but with an infinite linear chain standing for the bath.V) Two sites coupled both to different 
infinite chains (private environments). VI) Two sites coupled to the same infinite chain (public environment). 



Choosing the site energies with identical values: E^^n — ensures a continuum spectrum, while the election Ea = 
Eb — will further simplify the analysis. Two different alternatives for the spectrum dynamics arise when = 1 
(semi-infinite linear chain) and rij^ = —00 (infinite linear chain). The system-bath Hamiltonian Hg_£ depends on 
how we couple our two-site system to the environment (linear chain), but in general it will be of the form: 



Hs-s = ~Vo {c\c^i + cliC^ + c^gd^j + cljCgj , (4) 

where and fij label sites i and j in the environments v and fi respectively. 

It is crucial to stress that Vab , Vo a-nd V determine the relevant time scales of the whole problem. The first two give 
the rate of hopping from site A to site B and to the environment, respectively. The third is the jumping rate between 
sites in the environment. An "irreversible" decay to the environment, and hence the Fermi Golden Rule, implies 
that the unperturbed isolated system state has zero overlap with eigenstates of S + £. For this perturbation theory 
break down, the interaction with each environment eigenstate, Vo/-\/n must be much greater than the spacing between 
adjacent levels, of about V/N, i.e. the interaction time scale fi/Vb must be lower than the environment's Heisenberg 
time M/V. Imposing Vq <C Vab for all the cases we treated (weak coupling regime), we ensure the smooth degradation 
of the system's coherent evolution. When V ~ Vab the memory effects characterizing a non-Markovian situation 
lead to a very rich dynamical behavior. In the opposite limit, ii V ^ Vab then the validity of the FGR (Markovian 
situation) is expected to be recovered. Typically, this last situation will be represented in this work by a hopping 
V = 5Vab. 



III. NUMERICAL AND ANALYTICAL TOOLS 



A. SURVIVAL PROBABILITY AND LOSCHMIDT ECHO 



Two kind of measures for dynamical degradation are employed in this work: the Survival Probability (SP) PAA{t) 
and the local Loschmidt Echo. The SP is defined as 



PA^A{t)^ {A\exp ^iHt/h \A)Q{t) 



(5) 
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where 8(i) is the Heaviside step function, and the Hamiltonian H is defined by Eq. [T] It measures the probabihty 
of finding a particle in site A at time provided that the system has had a particle in the same site at time t — 0. 
In spin systems, this is a spin autocorrelation function (see Eq. IA3I and IA6I in Appendix |A| . The whole evolution of 
the system as reflected in the SP, is affected by a decay process, which is not trivial to separate from the intrinsic 
dynamics. Thus, to quantify decoherence, one relies on the observation of specific features as natural recurrences (Rabi 
oscillations or mesoscopic echoes) that appear at specific times. This limits the used time windows and limits the 
detailed assessment oi S — £ dynamics. Although it is not the perfect tool to quantify the effects of the environment, 
the SP behaves as a probe that reflects the overall dynamical process. 

With the purpose of get a continuous access to the S ^ £ dynamics that better reflects the environmental memory 
effects, we focus our attention on the LE. This measure has been used in the last years in different physical scenarios 
(both experimental [36r,39j and theoretical [s^ HO, i62i] ) in order to explain the behavior of the decoherence character- 
istic time. In general, the LE provides a direct measure of the decoherence process due to the environment, and even 
though it depends on the nature of the system's intrinsic dynamics, it does not depend much on its details. Its usual 
dynamical behavior presents an exponential decay regime [63j . which will be used to characterize the destruction of 
the system's coherent dynamics. 

The LE relies on the time reversal of the system's evolution and, in the present scenario, it has a direct physical 
interpretation. The LE can be understood as the measure of the amount of polarization returned to local site where 
it started. The controlled quantum dynamics is separated in two stages. First, the initial local excitation (particle 
in site A) evolves during a time ti, and then a time reversal procedure is applied during a period time {T — ti) 
which reverses the system's dynamics {Hg —Hg)- It is important to note that the bath's dynamics and the S — £ 
interactions {Hg and Hg-g respectively) remain unreversed during the backward evolution. This partial control 
results in a non-reversed perturbation E = {Hg + Hg-e) acting in both periods. Finally, the probability of finding 
the particle in site A forms a Loschmidt echo provided that ti = T/2: 



Mle{T)= (^|exp -i{-Hg + t){T -h)lh exp -i{Hg + t)h/h \A) 



(6) 



Clearly, in the case where the system is isolated, the local LE will have a steady value of 1. This means that the system 
is fully reversible. On the other hand, if the system is coupled to the infinite and continuous environment's spectrum, 
both the forward and backward effective Hamiltonians, Hg -\- S and —Hg + E respectively, become non-Hermitian. 
Thus, the LE decays with a characteristic rate, i.e., our reversal procedure has not been able to recover the excitation 
spread to the environment. The net decay of the norm of the state simply means the decay of the coherent part and 
then describes the non-trivial part of the Loschmidt Echo. 

For both survival probability and Loschmidt Echo evaluation, we diagonalize the Hamiltonian to obtain the evolution 
for every time. We use sufficiently large chains to approximate the nature of infinite ones. Indeed, the evolution times 
considered in this work are short enough to ensure the absence of dynamical finite-size effects (e.g. mesoscopic echoes 
appearing at the environment's Heisenberg time). A typical system is presented in Fig. [2]-a (the same model was 
deeply analyzed in Ref. f43|). 

We are interested in studying the decay rates as a function of the S — £ coupling parameter Vq (from Hamiltonian 
of Eq. HI). We fit the dynamics of the local LE with exponential decay function, which for the case plotted in Fig. 
[21-b, turns out to be the envelope of the Rabi oscillations in the SP and has the advantage of having a monotonous 
behavior. However, one should be aware that for an arbitrary system, the SP envelope does not necessarily match 
the local LE (the characteristic rates can be different). Every rate we obtain, is related to a particular choice of Vbj 
and we want to use them to asses the conditions of validity of a FGR regime: 



^^'^-l(Hg_y N,^ (7) 

n \ J 



where (^Hg^gj is a characteristic second moment of the S — £ interaction and A^i represents an appropriate density 

of directly connected states. Thus, we plot the decay rates as a function of Vq /hV (Fig. [2]-c). It is interesting to note 
that because of the linear chain topology this second moment coincides with Vq and A^i can be identified with a local 
density of state at the first site of the chain. In a general environment this correspondence can be assigned through 
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FIG. 2. (Color Online) a) Typical tight binding model considered, b) SP (dotted line) and the local LE (solid line) plots vs. 
time (in units of H/Vab)- c) Characteristic decay rate (exponential decay regime) as a function of Vq /hV (both in units of 
Vab/K). The points were obtained by the variation of the coupling constant Vo- 



B. ENVIRONMENT'S EFFECT IN A SELF-CONSISTENT FERMI GOLDEN RULE 



For the study of the analytical solutions we use the Green's Function (GF) formalism. In this framework it is 
possible to obtain the overall dynamics of the system by looking at the behavior of the GF poles [H, . Any element 
of the retarded GF can be obtained from the Fourier transform of the full propagator: 



oo 

lirn J {A\exp -i (^H - irji^ t/h \A) e{t) exp[+iet/h]dt 



s-Ea- S(e) e-Ea- A(e) + ir(£) ' 



(8) 



where E(£) is the appropriate self-energy with real and imaginary parts A(e) and — r(e) respectively. The self-energy 
operator is diagrammatically presented in Fig. [31 The bath memory is contained in their dependence on e, which arises 

from the bath's exact Green's Function Gii{e) at the directly connected site. Indeed, E (e) is given by |Vo|^ (e), 
see Fig. [3]-b. All the cases considered in this work (see Fig. [1]) can be reduced to the self-energy of a semi infinite 
linear chain [rii, = 1), i.e. 



S(e) = 



1/2 

iSL 

y2 



for |e| <2\V\. 



(9) 



For the infinite linear case (n^ — — oo), the previous expressions must be multiplied by 2. 

Notice that S(e) plays the role of the influence functional in the Feynman path integral formulation usually used 
to deal with memory effects for bosonic baths[63|. However, in such cases there are many more free parameters than 
in our case, i.e. the bath spectral density, the coupling strength with each mode and the temperature that fixes the 
bath occupation. In our spin environment model however, the last is simplified by the high temperature limit, while 
the two first become naturally determined by specific sum rules arising from the physical Hamiltonian we select. 

We remark that the complex GF poles are consequence of the unbounded nature of the system which prevents 
mesoscopic echoes [s^ and Poincare's recurrences. Also it is important to point out that we are working with a 
Hamiltonian in a parametric region holding a continuous spectrum where no localized modes appear. 
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FIG. 3. a) Diagrammatic representation for the retarded GF at site A, in the form of a Dyson's equation. The interaction 
with the environment is to infinite order in the self-energy given in (b) (see Eq. [8|. Simple lines with arrows are exact GF 
in absence of S — £ interactions, b) Self-energy diagram sums all orders in the hopping to the environment. Thick line with 
arrows is the exact propagator at a point of the isolated bath, denoted by Gn (see text). 

To get a better understanding of the connection between the dynamical behavior and the poles of the GF, we span 
the initial condition in the Survival Probability (see Eq. [S]) in terms of the energy eigenstates: 



PAAit) 



k=l 
oo 



eit) / de 



.fe=i 



exp [~iekt/h] 



The term in brackets is identified as the Local Density of States (LDoS) Nji(e) at the site A: 



(10) 



OC 

NA{e) ^ --lin J dtG'2^(t)cxp[-irf/?i] 



= -ilmG5^(£). 



(11) 



Therefore, we can identify Eq. [5] as the LDoS Fourier transform: 



PaA*) = 



OO 

e(t) J ^^Aie) exp[-iet/h] 



(12) 



The last expression can be numerically and analytically computed once we know the GF in the energy representation. 
As a matter of fact, since we have S(£), and hence Na{s), explicitly, we can compute the dynamics from the Fourier 
Transform mentioned. But, we prefered to obtain the numerical evolution (by exact diagonalization, as explained in 
previous section) of a finite environment. The interest in this last method arises from the fact that its use can be 
directly generalized to more complex systems (many-body), where an exact analytical solution is not accessible. The 
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analytical alternative based on the evaluation of the GF poles will enable to summarize the decay rate in a simple 
expression, as it is shown below. 

The characteristic decay of Pa.a(0 is determined by the bath's LDoS, iVi (e). This last can be obtained from the 

isolated bath's exact GF at site 1 (i.e. cfj^), 

N,{e) = --linGf,{e). (13) 

TT 

Hence, we recall that iVi (e) plays the most relevant role in S — £ dynamics as was quoted in Ref. [42]. In some 
systems the LDoS can be factorized as Na{£) = -^i(£) x Li{e) x L2(e), where Li{e) are Lorentzian functions (LFs) 
related to the real and imaginary parts (denoted by Ao and Fq respectively) of the GF poles. Hence, the convolution 
theorem applied to Eq. [pleads to a characteristic decay of Pa, Ait) ruled by Fq. In fact, the decay parameter of the 
exponential regime is given by 1/t — 2Fo/S. This is what we call self-consistent Fermi Golden Rule (SC-FGR) ^2§\. 

With the purpose of finding Fq, we focus on the Hamiltonian[Tl and follow the continued- fraction procedure described 
in Ref. For each type of S — £ coupling (see Fig. [1]), it is necessary to recalculate the poles of the GF. Due to 

the small number of poles of the systems, it was feasible to obtain the analytical solutions for all the cases presented 
in this work. Also, we address the behavior of Fq as a function of Vq <C 1. Indeed, we proceed to expand the solution 
near Vg — fl- 
it is important to notice that in the Taylor expansion, the linear and zero order terms vanish. Therefore, the 
imaginary part has Vq as the first non trivial term. As expected for a FGR, this is in strong agreement with Eq. [3 

2 

where we identify the order Vq as the second moment of the S — £ interaction (in general denoted by Hs-s )• In 

the next section, we present the corresponding values for 1/t = 2Fo/?i, expressed in the first non trivial order (as we 
said, the 2"^^), for each case considered. 

The usual alternative (easier and cheaper) to the presented scheme (SC-FGR) is the simple FGR, which is equivalent 
to evaluate the Green's function in a first pole approximation: 

[G^Aie)r' ^e^EA-mA) (14) 
^e-Ea- MEa) + iTiEA) (15) 



Since this yields an e- independent rate, it can be understood as a Wide Band Approximation (WBA) [66|-|68j. This 
approximation would imply neglecting any signature of dynamics and memory effects in the environment. Also, it 
misses some striking dynamical behaviors appearing at long times as the survival collapse [28, 42] and the subsequent 



power law decay [69|-|71) 



As a matter of fact, the WBA is represented by the condition of ^ Vab- Under this assumption, the environment 
acquires fast dynamics and the system does not receive any return from it (Markovian limit). In general, the common 
FGR has the form expressed in Eq. [T] There, the last factor (iVi) stands for the LDoS of the directly connected 
states, and by the application of the WBA, it is evaluated in the middle of the band spectrum (in our case, e = 0). 

In this work it is important to remember that the LDoS varies from the semi-infinite chain (rii, = 1, surface state) 
to the infinite chain (n^ = —oo, bulk state) as follows. 



A^i^.(£) = ^(^^ -^j e[2F -]£]]. (17) 

From the previous spectral structures, we stress the presence of van Hove singularities, which play an important role 
for long time dynamical behavior [l^li^, |6^. Indeed, the critical exponent characterizing the van Hove singularity 
can be related to the dimensionality of the space where the quantum excitation diffuses. Also, the two different 
convexities in the LDoS, near the middle of the band spectrum, will be of great relevance in the following discussions. 

IV. DECAY RATES: FGR AND BEYOND 

In this section we expose the main results of our work, summarized in Table |T1 The decay rates are presented 
as function of Vq/HV, for every case analyzed (Fig. [Ij and each approach employed (SP degradation, local LE 
degradation, WBA-FGR, and SC-FGR). 
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TABLE I. Degradation rates for all the methods analyzed (Survival Probability and local Loschmidt Echo degradation, Wide 
Band Approximation, and Self Consistent FGR). The cases are described in Fig. [T] If the system's and bath's time scales are 
equal, the physical situation is strickly non-Markovian. We point to restore Markovianity in cases where both scales differ by 
a factor of five. 



System 


SP degradation rate 


Local LE degradation rate 


WBA 




I 


2.04 ±0.05 


2.04 ± 0.05 


2 


2 


II 


1.00 ±0.02 


1.00 ±0.02 


1 


1 


III - {V = Vab) 


0.88 ±0.05 


0.88 ± 0.05 


1 


0.87 


III - (V = Wab) 


1.00 ±0.02 


1.00 ±0.02 


1 


0.995 


IV - = Vab) 


0.56 ±0.02 


0.56 ± 0.02 


0.5 


0.577 


IV - (1/ = bVAB) 


0.50 ±0.02 


0.50 ± 0.02 


0.5 


0.502 


V 


1.16 ±0.03 


1.16 ±0.03 


1 


1.15 


VI - (V = Vab) 


1.71 ±0.04 


1.20 ±0.04 


1 


1. 732 (/orward) and 0.577 (backward) 


Y1-(V = Wab) 


1.11 ±0.03 


1.02 ±0.03 


1 


I.106{forward) and 0.904 {backward) 



The rates predicted in the WBA column of Table HI correspond to the direct evaluation of Eq. [7l i.e.: 

l^^V,'N,,ie = 0), (18) 

where A stands for s (surface, for semi infinite chains) or b (bulk, in the case of infinite chains), in accordance to 
Eq. [inland Eq. [T7] respectively. Also, the details on the analytic calculation of the rates (SC FGR column, obtained 
by means of the GF poles) are presented in the Appendix |B1 

Systems I and II involves only one site. The difference between them is that I has a semi infinite chain acting as 
environment and II has an infinite one. The decay rates for these cases can be directly evaluated within the WBA, 
and agree exactly with the numerical solutions. 

From the cases III to VI, the system acquires its own time scale {H/Vab)- For these cases, we consider two time 
scales for the environment, as compared to the system's time scale. The first one, in which the time scales are equal 
{V = Vab) and the second in which the environment is "accelerated" by = 5Vab- 

Quite obviously, when the environment has the same time scale as the system, the WBA rate does not works at 
all. However, the rates computed by SP and LE degradation agree with the SC-FGR perfectly (except for the public 
environment, case VI, which will be analyzed below). This last quantity can be interpreted as the LDoS of the bath 
being evaluated in the exact solution for the eigen-energies of the whole {S and £) , and not in the middle of the band 
spectrum (see Fig. Hl-a). Indeed, the rate of the SC-FGR is obtained by evaluating the LDoS at the real part of the 
GF poles. 

It is important to notice that the rate obtained for the SC-FGR are higher or lower than the predicted by the 
WBA, depending if the bath is infinite or semi-infinite. We can interpret this result by looking at Fig. 21 Once we 
evaluate the LDoS in Re (Spo/e) — ^Vab/V, we observe that the values are lower for the SC-FGR than the WBA in 
the semi-infinite LDoS, otherwise higher if the infinite LDOS is considered. 

At this point, for the two level system, we are able to link the convexity of the bath's LDoS to the observed decay 
rate. We notice that if the LDoS is convex, the exact rates are greater than the WBA prediction, but instead if it is 
concave, the rates are smaller. In other words, it depends on the shift away from the middle of the band spectrum, 
moving along the LDoS slope. 

In order to point to the strictly Markovian case, where any return from £ to 5 is suppressed, we examine accelerated 
environments. For these cases, all the rates converge to the same value. From Fig. |3]-a we can see the convergence of 
the exact solutions towards the middle of the band spectrum, as long as the condition V Vab is better fulfilled. 

So far we have only considered a single and private environment, and in the following cases we address the private- 
public discussion. As we mentioned before, private baths act over each site individually, and a single public bath acts 
over two (or eventually more sites) of the same system. This public S — £ interaction induce new types of correlations, 
increasing the dynamical complexity of the physical process. 

For case V we explore the possibility of two private baths (i.e. two infinite environments each one connected to 
one site). Since here both sites, instead of only one, are affected by the environment, it is expected a double rate 
compared to the case IV. In Table U we found that the simulation results agree with the analytical prediction and 
again, the rates are greater than the WBA value, which means that we are moving with the exact eigen-energies 
through a convex LDoS, away from middle of the band. 

Finally we analyze a very interesting case, where the system is in presence of a public bath (VI). For these cases it 
is observed that the SP rate differs from the local LE. This behavior shows certain asymmetry between the forward 




FIG. 4. (Color Online) Local Density of States as function of energy (in units of 1/Vab and Vab respectively) for: a) Semi- 
Infinite chain (Surface) and b) Infinite chain (Bulk). The horizontal dotted lines are the values Vab/V which stand for the 
system energies (renormalized with V). The vertical dotted lines together with their numerical values are the LDoS for each 
system in the normal {V = Vab) or the fast {V — 5Vab) configuration. 



and backward evolutions. Using the SC-FGR approach it can be shown that the imaginary part of the poles depends 
on the relative value of Vab and V (see Appendix [B)) . 



i_2 ^^ yw^-v^^ v,' 

This dependence induces a different decay rate whenever the system evolves forward {Vab > 0) or backward {Vab < 0). 
Thus the SP and the local LE are not equivalent (see Table It is worth mentioning that the asymmetry (dependence 
on the relative sign) in the rates for the forward and backward evolutions, arises only when the bath is public. 

The analytical prediction for the local LE rate is obtained observing that the total evolution for reversed dynamics 
is proportional to a product of two exponential evolutions (forward and backward), 



Mle{T = 2t) cx exp(- — ) exp(--) 

Tf n 



exp 



TfU 



exp 



-T 



V + n 



(20) 
(21) 



where l/r/ and 1/ti, correspond to the forward and backward degradation rates (see also case VI in Table |T|. Eq. EOl 
shows that the mean rate for a single time-reversed evolution should agree with the LE decay {1/tle = 1-20 ± 0.04 
VQ/hV). In fact, this is true, since the two SC-FGR rates (l/r/ and l/n) yield a mean value 1/t = 1.15 Vq/HV. 
Moreover, and quite obviously, the forward rate agrees with the SP decay. 

Let us now link the decay process to an appropriate LDoS. As it is clearly shown in Appendix [B] we can transform 
the original model of case VI to an equivalent one (see Fig. [5]). If a suitable change of basis is applied (basically 
turning into symmetric and anti symmetric basis), then it is only necessary to analyze two semi infinite linear chains 
(a particular case treated in Ref. 0, [llj). Since the initial condition has equal weight on both effective chains, the 
corresponding rates for them have to be added. Further details on the symmetrization transformation are explained 
in Appendix [Bl 

From the LDoS considered in Fig. [6] we can explain why the rate l/r/ is always above the WBA limit, and the 
l/rf, rate is always below it. Moreover, as long as the forward rate agrees with the SP decay, this also supports the 
SP rate being greater than the WBA value. 

Once more, we make the bath's dynamics faster, leading to a better applicability of the WBA and the rates converge 
again to the traditional FGR description (middle of the band spectrum). 

Now we briefly analyze the characteristic local LE decay times for a "highly public" bath (See Fig. [7]). This case 
is of interest for studying spin dynamics in ring and ladder like systems (see Ref. (32j). 
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FIG. 5. (Color online). The schematic symmetrization procedure of the 2-sites system interacting with the public environment. 
The original problem can be cast as two independent semi infinite linear chains. 
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FIG. 6. (Color online) LDoS as function of energy (in units of 1/Vab and Vab respectively) for the equivalent problem of the 
public case, after symmetrization transformation. Forward and backward stages of the evolution are indicated. The horizontal 
dotted lines are the values Vab/V which stand for the system energies (renormalized to V). The vertical dotted lines together 
with their numerical values are the LDoS for the system VI in the normal {V — Vab) or the fast {V = 5Vab) configuration. 
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FIG. 7. (Color Online). Fermion model. A five-site system laterally coupled to an infinite chain. 



In this system we will consider three different relation for time scales: V — Vg, V — 5Vs and V = 8.75T4 (we also 
assume V Vq). For these cases the decay rates (local LE degradation rate measured in units of Vq/HV) computed 
by the numerical solution are (2.66 ± 0.07), (1.54 ± 0.05) and (1.16 ± 0.04), respectively. Therefore, we observe that 
the first two results do not match with that expected by WBA. Moreover, as we have seen before, the convexity of 
the bulk LDoS, induce bigger values. By setting V — 8.7514 (approximately one order of magnitude of difference) we 
are near the WBA validity and thus nearly reproduce the expected result: an energy independent rate, of value 1.00 
Vq /HV, in accordance to Eq. [181 We stress here that the strong sensibility on the relation between time scales, is 
enhanced by a highly public S — £ interaction. 



V. FURTHER DISCUSSIONS AND CONCLUSIONS 



Traditionally, the first approach to evaluate decay rates is the FGR based on an underlying WBA. We have shown 
several variations of a simple SWAP gate interacting with an environment (or several ones) where this scheme is not 
as good as one might assume. This failure, which can be understood as non-Markovianity, is due to two contributions 
(which should not be considered as well separated effects). 

First, if the WBA limit {V ^ Vab) is not well fulfilled, then it is quite obvious that the energy independent rate will 
not be representative as it would oversimplify the decay process. When the system's time scale becomes comparable 
to that of the environment (i.e. Vab — V, non-Markovian situation), the decay rate departs from the usual FGR 
approximation evaluated at single energy level. Depending on the bath's spectral structure (LDoS) the actual rate 
can be greater or lower. We have shown that if the LDoS is a convex function of the energy, then any shift away from 
the middle of the band spectrum will produce a greater decay rate. But if the LDoS is a concave energy function, the 
rate would be lower. In general, we have seen that dynamical complexity (understood as how much the decay departs 
from the WBA validity) arises when the time scales of the system, the environment and the S — £ interactions, are 
commensurable. 

Second, when the dynamics is non-Markovian due to similar time scales, the possibility of interaction through 
the environment becomes appreciable. In a public bath, the correlations generated by the multiple connections 
between different parts of the system and the bath, are more effective to depart the physical process away from the 
WBA validity, for a given relation of time scales. In general this leads to too subtle correlations to be treated in 
spectral models for the bath which here becomes c^uite natural through the use of specific Hamiltonian models of the 
environment. 

In order to address the question on the private and public nature of the S — £ interaction, we confront the cases of 
two independent baths (case V) and one common bath (case VI). In cases where S — £ interaction is private there is 
no dependence of the rate on the relative sign of V and Vab- On the other hand, if the S — £ interaction is public, 
the forward and backward evolution will change the decay rates and will produce a mean value for the local LE 
decoherence rate. It is important to stress here that we are comparing cases where S and £ have similar time scales. 
In this condition, it is observed that decay rates for the public case are greater than for the private. In this sense 
we interpret (at least for this cases) that the private bath is less harmful than a public one. Strictly speaking, this 
is a consequence of moving along the LDoS in the effective FGR (see Figs. |4]and[6]). A public scenario produces an 
enhancement in the alteration of the FGR away from the WBA limit, as compared to the private one. As a matter of 
fact, private bath is less efficient for correlated memory-like returns to the system. As it has been previously stated in 
a general sense, and based on complete formal grounds, the feedback of information from £ to 5 is the central issue 
for quantifying a non-Markovian bath [t^I • 

On the other hand, we have shown how Markovianity is restored by changing the relation between S and £ time 
scales. In fact, for all the cases treated in this work, we have found that once we move towards the WBA limit 
(represented by a bath dynamics 5 times faster than the system's), the effects of a public or private bath are no longer 
relevant. For both cases the environment is fast enough to wash out all the memory effects, since any inner excitation 
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is rapidly spread. At this point the rate for both situations become equal and agree with the one predicted by the 
common FGR. This means that when the system's and the bath's time scales are well differentiated (specifically when 
the bath's is very short as compared to the system's), then the public-private reservoir discussion becomes irrelevant. 

In many physical situations, when addressing 1-D and l-D"*" systems (see, for example [s?]), a coupled environment 
behaves with a convex LDoS. Therefore, memory effects and complexity in the structure of the S — £ interaction 
produce an enhancement of the coherent dynamics degradation. Additionally, we stress here the wider "spectral 
exploration" of the Loschmidt Echo as compared to the Survival Probability (a difference clearly shown in the public 
S ~ £ interaction case, where the global LE rate is a mean value of two non symmetric processes). 

Even though we do not claim full generality for the results discussed here, there are some interesting universal issues 
to remark. The use of the Wide Band Approximation (simple FGR) can easily lead to quantitative and qualitative 
errors if the time scales and the physical structure of S and £ are similar. Furthermore, the way they are coupled 
to each other [public or private), plays a fundamental role. In general, dynamical complexity grows when those 
characteristic times are similar and when there is no privacy in the S — £ interaction. 
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Appendix A: Excitation dynamics in 1-d systems and the Spin-Fermion mapping 



In previous sections we developed a scheme that could be used to treat spin polarization dynamics [3l|, |33| , under 
certain assumptions. In that sense, the well-known mapping between spins and fermions [sij ] has been used for 
formulating the spin problem in terms of the non-equilibrium Keldysh formalism [sol Is^ . We briefly present here 
the Spin-Fermion mapping, and how the tight binding models discussed along the present article arise. 

A simple case that can be treated in this context is a linear chain of M spins in an external magnetic field. They 
interact with their nearest neighbors through XY coupling: 



M-2 



M-1 M-2 ^ 



(Al) 
(A2) 



Tl = 



ri=0 



where are the rising and lowering operators S*^ = S*^ ± iS^. The dynamics of the Af-spin system, evolving 
under the Hamiltonian H, is usually described by means of the two site spin correlation function. 



.,|^|(to)Sf(to)|*e,) 



(A3) 



The quantity of Eq. I A3 1 gives the amount of local polarization in the z component at time t on the /th site, provided 
that the system was, at time to, in its equilibrium state with a spin up added at ith site. Also, S^f{t) = e'^'S'fe 



{N)\ 



-iHt 



eq 



I 



is the thermodynamical many-body 



is the spin operator in the Heisenberg representation and \^eq) — '^n 

equilibrium state constructed by adding states with different number N of spins up with appropriate statistical weights 
and random phases. 

The Jordan- Wigner transformation (JWT) links spin and fermion operators at each site, by the following relation: 



exp 



i^^ 



At p 
mm 



m— 1 



(A4) 



where c\^, Cm are the canonical fermionic operators. The use of the JWT on the Hamiltonian IA2] vields: 
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M-l , A/-2 



n=0 n=0 



where e„ = hftn are the site energies and Vn+i,n = ^Jn+i,n are the hoppings. Due to the short range interaction 
(first neighbors), after the apphcation of the JWT, the only non-zero coupling terms between spins are proportional 

to cJj_(_jC„ = S^^-^S~. Each subspace with (^) states of spin projection ^nj = N — M/2 is now a subspace 

with N non-interacting fermions. The eigenfunctions 'I'^^'^ are expressed as a single Slater determinant built up 

with the single particle wave functions (pa of energy Under these circumstances, and setting |i) = c\ |0) (with |0) 
the fermion vacuum), Eq. IA3I reduces to: 



Pf,,(t) = (/|exp[-ii?i/ftj \i)Q{t) 

where G^^{t) is the retarded Green's function for a single fermion that connects sites / and i. 
While similar steps lead to description of excitations in Double Quantum (DQ) Hamiltonian [5^ 



(A6) 



Hdq — 



M-2 

E 



J n 



(A7) 



those will not be detailed here. The fundamental issue is the underl yiiig on e- body dynamics, which for the DQ 
Hamiltonian is revealed by a unitary transformation H^q = U^HxyU 13, 5^ that links it to an XY Hamiltonian. 

Hence, for one dimensional (1-D) chains of spins with first neighbors XY or DQ interactions, and in the high tem- 
perature regime, the dynamics of an excitation (either an injected local polarization or multiple quantum coherence) 
is completely equivalent to the evolution of a single particle wave function, ruled by a tight-binding Hamiltonian. 
Therefore, it turns out that the analysis of Hamiltonians like the one in Eq. IA5[ can be casted for treating and 
studying several effects in spin chains [ill (a typical scenario in QIP). 



Appendix B: Green's Function Poles. 

Here we present the detailed analytical derivation of the GF poles and further approximations. The results from 
this appendix have been summarized in the Table U showed in Sec. IIVI 

The first model under consideration (Fig. [T]-I) is given by one site coupled to a semi-infinite chain. In this case the 
GF pole results purely imaginary. 

Even though this is indeed the exact solution, we span it for Vq /V <C 1 and obtain: 

Epoie — ■ (B2) 



Thus the theoretical decay rate for this model can be expressed as: 



i4ro4M.,..)4f, ,B3) 



If we change the environment to an infinite chain (Fig. [TJ-II) the pole is. 
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1/2 



^poU = -i[2VHl-^§, + l)j --if (B4) 
The case III, which corresponds to two sites coupled to a semi-infinite chain, has been solved previously in Ref. p^ . 



Using the solutions of Eq. (26) and (27) in Ref. we can directly evaluate the real and imaginary parts of the 

poles, 



4 



1/2 



±Vab[1-^^] (B6) 



W V 



Tn = 



1/2 



2 



1 y/w^ ~ vis 

^ V V 



(B7) 



Now, we change slightly the geometry of these systems and consider again an infinite chain as environment instead 
the semi-infinite. For these cases we will write only the first non trivial terms of their Taylor expansion instead of the 
full solution. For the model of Fig. [UIV, the solution is: 

Ao ^ ±Vab + O(^) (B8) 
Fq--- - — 1 — - ^-^^nrJL\ (B9) 

The next case (Fig. [T]-V) is a two-site system coupled to two private baths. In this case, the solution is, 

Ao~±VSis+0(^) (BIO) 
V V^ v^ 

Tr+0(T7J)- (Bll) 



^W^-Vls V 

At this point we observe that the imaginary part for the case V is twice of the system IV. This behavior is consistent 
with the count of the "number" of private baths connected to the system, and the proportion affected by those private 
baths. 

Finally, for the case of Fig. [1]-VI, the model involves a two-site system coupled to a common bath. For this case 
the solution is expressed in the following form, 

^^^±\Vab---^)+0{^) (B12) 



^° - " W-2VAB V" + ^^V^^- ^^^^^ 



Notice that the imaginary part of the pole is linear in Vab in a sign dependent manner (in contrast to previous 
cases). This linearity translate into a different value for Fq depending on the sign of Vab- As a matter of fact, the 
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local LE indeed relies on the change of sign to revert the dynamics. To understand the physics below this difference 
we have to identify the LDoS involved in the decay process. Accordingly, we symmetrize the basis, thus dimmerizing 
the system, as shown in Fig. [5] (i.e. we take pairs of site states and map them into symmetric and anti symmetric 
states.) Therefore the tight binding Hamiltonian for this case is mapped to the form of Eq. IB14I 

l^> \B) |1> 1-1) |2) 1-2) 

{A\ ( Vab Vb ^ 



{B\ 

(1| 

(-11 

(2| 
(-21 



Vab 

Vq V 
V 



V V 



V 



V 



J 

\2s) lis) \ABs) \ABa) \Ia) |2a) 



(2s| 
(Isl 

{ABa\ 
{U\ 



V 

V V Vf) 
K, Vab 



-Vab Vo 
Vo -V V 

V 



(B14) 



where \ns) = (|n) + j-n)) /A |n^) = (|n) - |-n)) /V2, n = l,2,3..and \ABs) = {\A) + \B)) / \ABa) = 
{\A) — \B)) / a/2. From the Hamiltonian H' (Eq. IB14[) it is easy to identify the splitting of the original problem into 
two semi infinite tight bindin g ch ains, with only the first two site energies non zero. This problem has been previously 
addressed in Refs. and (72l| and we recall the LDoS computed there. Thus in Fig. [Blwe identify the energy 
excitation on the spectral structure of the environment, which is relevant during the decay process. 
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